Anomalous growth of thermoelectric power in gapped graphene 
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There exist experiments indicating tliat at certain conditions like an appropriate substrate a gap 
of the order of 10 meV can be opened at the Dirac points of a quasi-particle spectrum of graphene. 
We demonstrate that opening of such a gap can result in appearance of a fingerprint bump of the 
Seebeck signal when the chemical potential approaches the gap edge. A magnitude of the bump can 
be up to one order higher than already large value of the thermopower occurring in graphene. Such 
a giant effect, accompanied by the non- monotonous dependence on the chemical potential, is related 
to the emergence of a new channel of quasi-particle scattering from impurities with the relaxation 
time strongly dependent on the energy. We analyze behavior of conductivity and thermopower in 
such system accounting for quasi-particle scattering from impurities with the model potential in self- 
consistent scheme. Reproducing the existing results for the case of gapless graphene we demonstrate 
a failure of the simple Mott formula in the case under consideration. 

PACS numbers: 65.80.Ck, 81.05.ue 



Control of heat flows and minimization of heat losses 
is the important aspect of designing of modern nano- 
electronic devices, in particular, based on graphene 
The experiments indicate Q that thermoelectric effect in 
graphene accounts for up to one-third of the contact tem- 
perature changes and thus it can play a significant role 
in cooling down of such systems. The measured ther- 
mopower S reaches the values ks/e IQO fiV/K at room 
temperatures, where fc^ is the Boltzmann constant and 
— e < is the electron charge. This and other experimen- 
tal results on thermoelectric transport in graphene 
were understood theoretically mostly basing on a sim- 
ple Mott relation between thermopower and logarithmic 
derivative of the electrical conductivity a{fi, T), 

S{f,,T) = -^T^[\na{f,,T = 0)], (1) 

where ^ is the chemical potential, T is the temperature 
(we set ks = 1) 0- However, the cited experiments show 
that the Mott formula ([T]) fails when jj, approaches the 
vicinity of the Dirac point at high temperatures, espe- 
cially in high- mobility graphene [6|. The available theo- 
retical analysis (sl-fioj shows that the failure of Eq. ([T]) can 
be attributed to breaking of the conditions of its appli- 
cability, which read as T <C and/or T <C 7 (here 7 is 
the characteristic energy scale on which the conductivity 
a{fj,,T = 0) varies around the Fermi level). 

The purpose of the present letter is to show that 
already large value of the thermopower occurring in 
graphene can be further increased up to one order of 
magnitude by opening in this or that way (see e.g. 
Refs. [m [lH) of a gap A in its quasi-particle spectrum. 



We will show that such opening is accompanied by the 
emergence of the new channel of quasi-particle scattering 
with the relaxation time strongly dependent on energy 
and it results in appearance of the giant bump in ther- 
mopower when chemical potential approaches the gap 
edge. The situation here turns out to be very similar 
to the well known anomaly of thermopower close to the 
electronic topological transition (see Ref. [l^ 14 1 for a 



review) related to the scattering of electrons from all ex- 
tended periphery of Fermi surface to the "trap" , pre- 
sented by the small new void or narrow "neck" of the 
latter. 

It is worth mentioning that some experimental findings 
isl . 16] indicate that, indeed, a gap in graphene quasi- 
particle spectrum opens at the Dirac point, and probably 
it can be attributed to the effect of a substrate. Yet, 
this issue still remains an open problem of physics of 
graphene |lll.ll2l|. This is why, in principle, one can turn 
the obtained below results on the contrary and use the 
thermopower measurements as a sensitive probe of the 
gap opening in graphene spectrum. 

The Hamiltonian for graphene can be written down in 
the momentum representation as 



H = 



where 



E / 7|\n(p)[i^(p)-^?„]T.(p), 



n{p) = ?+cj) (p) + (p) + Afg, 



(2) 



To, T3,r± = (ti ± iT2)/2 are Pauli matrices acting in 
the sublatticc space on the spinors, To-(p) and Tj.(p) = 



2 



(a^(p), 6j.(p)) with the creation (annihilation) operators 
of electrons at(p)j ^tip) (QctIp): ^cr(p))j corresponding 
to A and B sublattices, respectively, and spin subscript a. 
The full form of the complex function 0(p) is provided, 
for example, in Ref. [llj and in present consideration it is 
important only that near two independent K points the 
dispersion ^ = |0(p)| = ^iifIpI, where vp is the Fermi 
velocity and the the wave vector p is measured from cor- 
responding K point. In result of diagonalization of the 
operator ([2]) one finds that presence of the gap A in it 
breaks the equivalence between A and B sublattices, and 
the spectrum of the quasi-particle excitations close to K 
points takes the form E (p) = zL^/Wvj/p^+A^ — fi. 

Quasi-particle scattering from impurities we will ac- 
count for in the framework of the Abrikosov-Gorkov 
scheme writing the self-consistent equation for self-energy 
in the matrix form 



S(p,e„) 



BZ 



^(q)G(p-q,£„)F(q), (3) 



with Hi as concentration of impurities and e„ = 7rT(2n-|- 
1), the full inverse Green's function 



G-i(p,£„) = Go'(p,e„)-S(p,e„ 



and the free inverse Green's function 



Go ^(P,£n) = {i£n + f-l)To - nip). 



(4) 



(5) 



The sublattices A and B are shifted for the distance of 
the order of lattice constant a, hence their images in in- 
verse space are separated by the momenta of the order 
of h/a. Hence, for relatively long-range potential V{q) in 
Eq. ([U one can ignore the quasi-particle scattering be- 
tween the inequivalent valleys. At the same time we will 
assume V{q) as momentum independent for the intra- 
valley scattering, i.e. 



^(q) 




I I ^ max{|/x|,A} 



max{ , A} 



< ^ 



(6) 



Next step in solution of Eq. ([3|) is the decomposi- 
tion of the self-energy over Pauli matrices E(p,e„) = 
^^^0 (Ji(p, e„)fi. One can see that ignoring self- 
consistence procedure in Eq. ([3]) with potential ([6]) (which 
does not mix valleys), one obtains S'°^(p,e„) in the di- 
agonal form. It is possible to show that the off-diagonal 
components appear in the order nf only, i.e. crf-(p,£„) 
and (T^(p,e„) terms can be omitted. Finally, the mo- 
mentum dependence of I](p,£„) can be also ignored in 
view of the constancy of potential (|6]) in the domain of 
each valley and taking into account that the main contri- 
bution to E appears from the relatively large momenta 
q| ~ max{\ij,\, A} /Hvp- This latter fact justifies our 
use of the Abrikosov-Gorkov technique for averaging of 
the quasiparticle scattering over impurity positions. In 



result, after the analytical continuation ie„ — >■ e, the ma- 



trix Eq. (jSj reduces to the system of two equations 17[: 
4h 



X 



w 



(7) 



/ [e + fi- a§ ie)]-e-[A + (e)] ^ 
where we introduce the "relaxation time scale" [l6 
l_ ^ n,\u{Q)\^\^J.\ 

To 



(8) 



The real parts Re and Re which logarithmically 
dependent on the high energy cut-off W, can be included 
to the renormalized fj, and A, respectively. For the self en- 
ergy imaginary parts, which determine the quasi-particle 
relaxation rate, one finds: 



ImcrQ'(£) 

h 
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Im crfie) 



A^ 



sgn(£ + pi) 



(to/^)|mI 



e + A* 
A 



(9) 



One can note that the denominator of the full GF (jlj 
can be approximated as: 

[e + n- ilmaff (£)]2 - - [A + ilmaf (e)]^ 
^{e + ti + iT{e)f -E\p), 

so that the full inverse GF can now be written as 
[G^'(p,e + iO)]-i « [G^'(p,e)]-i + i%T{e). Here the 
energy-dependent scattering rate r(e), central for our 
consideration, explicitly appears: 



T{e) = -Ima^ie) - 



£ + /i 

A2 



I^ + mIImI 



(10) 



with Fq = 2h/TQ. In the numerical results presented 
below we use the value Fq = 20 K ignoring its dependence 
on the carrier concentration. 

It follows from Eq. ((TUl) that for {e + < A^ the 
scattering is absent. Further consideration shows that 
in spite of the presence of F(e) in the denominators of 
Eq. (fT2)) . this fact does not result in divergency of the 
physical observables. Yet, one should keep in mind that 
some scattering processes beyond our model along with 
the next order corrections to the solution © can make 
the scattering rate finite below the gap edge. For our 
numerical work we took this into account by adding a 
small residual scattering rate 70 to r(e). In accordance 
to the theoretical analysis the final results turn out to be 
practically independent on the value 70. 

Using the Kubo relations one can derive electric con- 
ductivity and thermoelectric coefficient in the explicit 
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form 



deA{e,T{e),A) j 1 

e/{eT) 



2rcosh 



2 _^ 

2T 



(11) 



where in the presence of A the function A is given by 



^(£,r(£),A) 



1 

2^ 



(A^ + £)^-A^ + r^(£) 

2|^ + e|r(£) 



X I — - arctan —-— 

2 2|^ + e|r(e) 



(12) 



For A = Eq. (1121) reduces to the commonly used ex- 
pression (see e.g Refs. [1, [13] )■ In this case, setting 
also r(£) = Fq = const, one obtains for |//| ^ TjTq 
that a = e^\u\J(2'KhTo) and /3 = 7rersgn/i/(6SFo), in 



agreement to [21|. Then the value of the thermopower 
S = —13 /a turns out to be the same as in the conven- 
tional metals, S = — (7r^/3e)T///, and coincides with the 
result obtained directly from the Mott formula ([T]). 

The dependences I3{^i), and S'(/i) are shown in 

the top, middle, and bottom panels of Fig. [U respectively. 
The left side [Fig. [U (a)] is for T = 1 K and the right side 
[Fig. [U (b)] is for T = 5 K. The dashed (red) curves 
in all panels correspond to the reference case A = 0, 
r(e) = Fo, so that oc j^j and S{y) cc 1/fi for large 
1^1 . The general expressions pT|) and ([T^ allow also to 
reproduce the gapped case with the energy-independent 
scattering F(£) = Fq 0, 21 1. The corresponding refer- 
ence dependences computed for A = 50K are presented 
in all panels by the dash-dotted (black) curves. 

Our main results obtained with the energy dependent 
r(£) given by Eq. PU)) and A = 50K are presented in all 
panels by the the solid (blue) curves. 

The behavior of the conductivity a{fi) is shown in the 
top panel of Fig. [1] (a) and (b). One can see that it dras- 
tically changes due to account for the energy dependence 
of F(£). Namely, a clear kink in the dependence cr(fi) ap- 
pears at the gap edge, at |^| = A, while below it, when 
< A, the value of cr(/i) is strongly suppressed as com- 
pared to the case F(£) = Fq, A = const. This kink is 
smeared out with the growth of temperature [Fig. [T] (b) , 
T = 5K]. 

Let us pass to the discussion of the Seebeck signal. 
In the gaplcss case with the constant scattering rate 
F(£) = Fq the signal monotonously changes with fi pass- 
ing zero without any visible anomaly. Account for gap in 
the quasiparticle spectrum shows up itself by smoothed 
bumps at \fj,\ = A, which are rapidly smeared out by 
temperature [middle and bottom panels of Fig. [T] . 

The behavior of thermoelectric coefficient (3{fi) and 
thermopower S{fi) in the case under consideration of 
gapped graphenc with the energy-dependent relaxation 
time is shown in the middle and bottom panels of Fig. [1] 
by the solid curves. Let us stress that the curves corre- 
sponding to S{fi) in the bottom panel are multiplied by 
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FIG. 1: (Color online) Top panel: electrical conductivity a 
in units of the value ctq = 2e^ / {■k^K); middle panel: thermo- 
electric coefficient /3 in units of the value j3o = kse/h; bottom 
panel: thermopower S in units of the value ks/e as functions 
of the chemical potential ^. Left side: (a) for T = 1 K and 
right side: (b) for T = ^K. In all graphs three cases are 
shown: the dashed (red) curve is for the energy independent 
scattering F = To and A = 0, the dash-dotted (black) curve 
is for F = Fo and A = 50 K, the solid (blue) curve is for 
energy-dependent r(£) and A = 50K. The solid curve in the 
bottom panel is multiplied by the factor 0.2 The thin lines on 
the right bottom panel are obtained using the Mott formula. 



the factor 0.2 to present them conveniently on the back- 
ground of the previous cases. This means that the peaks 
of the Seebeck signal arc at least 5 times higher than 
those ones obtained for F(£) = const case. 

A strong enhancement of the Seebeck signal in the case 
when F(£) is energy dependent one could already foresee 
even basing on the Mott formula ([l}. However, this for- 
mula gives only a hint on the singular behavior of the 
Seebeck signal and cannot be used for any quantitative 
description. Indeed, the thin lines in the bottom panel of 
Fig. [T] (b) are computed using the zero temperature elec- 
trical conductivity a{^.,T = 0) = (2eV^)^(0, F(0), A) 
and the Mott formula ([T]), while the thick lines in the 
middle and bottom panels of Fig. [T] are plotted using the 
Kubo formulas ([TT|) both for a and /3. One finds that for 
the case of A = and F(£) — const agreement between 
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the Kubo and Mott formulas is very good. We checked 
that for T = 1 K it becomes perfect, so that the hnes for 
the Mott formula are not shown in the left bottom panel 
of Fig. [T] The right bottom panel shows that for the case 
of finite A and r(e) = const one can already see some 
discrepancies between the Kubo and Mott formulas, es- 
pecially near \fj,\ = A. Finally, the Mott formula fails 
completely when the energy dependence of r(£) is taken 
into account. 

Two more comments on the obtained bump of the See- 
beck signal should be done. First, the shape of the bump 
depends on the presence ~ A^ term in Eq. (fTO|) and 
thus account for the self-energy (73 (e) is important for 
the qualitative theory. Second, our arguments are also 
directly applicable to a gapped bilayer graphene, and in- 
deed the computations done in [22j confirm this. 

The applicability of the model potential © to the 
case under consideration deserves a more detailed dis- 
cussion. It is worth to stress that we use it to solve 
the equation for self-energy with further fixation of Fq = 
2hTo^{\li\=A) = n,\u{0)\^\A\/{2h^vl) [see Eq. 
This procedure gives the consistent analytical treatment 
of the problem close to the gap edge |^ — A| <^ A, but 
it does not allow to reproduce correctly the experimen- 
tally observed dependences of a and S on the carrier 
concentration n(oc /i^sgn/x) beyond this region. In or- 
der to get the better agreement with the experiment in a 
wider interval of concentrations one could use the scatter- 
ing potential V{q) in the form of a long-range Coulomb 

UM " 



one (see, for instance, Refs. [9|, [101, [IE 
consideration one obtains that at large n the scattering 
rate Fq cc (contrary to our reference case, where 

Fq = const), that results in the observed linear depen- 
dence cr(n) oc \n\ (contrary to our a{n) cc 

The specifics of thermopower consists of its sensitivity 
to derivative of the scattering rate. This is why presence 
of the step function in Eq. (|10p produces much stronger 
effect on the behavior of S{fi) in the vicinity of \^\ w A 
than a relatively slow energy dependence of Fq which 
could appear from the screened Coulomb potential V{q). 
Let us attract again the readers attention to the evi- 
dent analogy between the transport in gapped graphene 
and that one in metal being close to electronic topolog- 
ical transition. Indeed, in vicinity of the critical point 
fi = fic, when the Fermi surface connectivity changes, 
the quasiparticle relaxation rate also acquires the contri- 
bution depending on energy in the form of step function, 
what generates the well known kinks in conductivity and 
peaks in thermopower [l^ . 

One can imagine that designing future nano-electronic 
devices it will be possible to control their temperature 
regime using Peltier cooling effect which is also governed 
by the value of the thermoelectric coefficient. In this 
respect easily tunable band gap in a bilayer graphene 
looks more promising for practical applications. 

Turning this around, one can exploit the predicted gi- 



ant peak of the Seebeck signal as a signature of the gap 
opening. Its existence in the quasi-particle spectrum of 
single layer graphene presents an intersti ng pro blem. Al- 
though there are not so many evidences 15|, ll6[ that this 
gap is present in zero magnetic field, there is a growing 
confidence that = quantum Hall state in graphene is 



gapped (see Ref. [2J] and references therein), so that a 
generalization of the present work for a finite magnetic 
may present some interest. 
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